Mon. Not. R. Astron. Soc. OOP. [HO (2002') Printed 2 February 2008 (MN I*TeX style file v2.2) 



On the effective particle acceleration in the paraboloidal 
magnetic field 



V. S. Beskin^* and E. E. Nokhrina^ 



o 
o 

(N 
> 

o 
O 

(N 
> 

cn 
m 
\o 

o 
in 
o 

O 



^ P.N.Lebedev Physical Institute, Leninsky prosp., 53, Moscow, 119991, Russia 

^Moscow Institute of Physics and Technology, Dolgoprudny, Moscow region, 141700, Russia 



Accepted 1988 December 15. Received 1988 December 14; in original form 1988 October 11 



ABSTRACT 

The problem of the efficiency of particle acceleration for paraboloidal poloidal magnetic 
field is considered within the approach of steady axisymmetric MHD flow. For the 
large Michel magnetization parameter a it is possible to linearize the stream equation 
near the force-free solu tion and to so lve the problem self-consistently as was done by 
iBeskin. Kuznetsova fc Rafiko\^ l)1998f) for monopole magnetic field. It is shown that 
on the fast magnetosonic surface the particle Lorentz factor 7 does not exceed the 
standard value a^^^. On the other hand, in the supersonic region the Lorentz factor 
grows with the distance z from the equatorial plane as 7 « {z/Ri^y^^ up to the 
distance z « (t^R-l, where i?L — c/Op is the radius of the light cylinder. Thus, the 
maximal Lorentz factor is 7max ~ cr, which corresponds to almost the full conversion 
of the Poynting energy flux into the particle kinetic one. 
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1 INTRODUCTION 

An activity of many astrophysical sources (pulsars, active 
galactic nuclei) is associated witli the presence of a strong 
magnetic field (~ 10* G for AGNs and ~ 10^^ G for pulsars) 
surrounding the rapidly rotating object and a relativistic 
particle outflow. Convenient way to characterize such flows 
is to int roduce t he ma gnetization parameter a, as was first 
done bv lMicheJ (Ec 



(1) 



Here ^Ptot is the total magnetic flux, and A = n/nGj is the 
multiplication parameter of plasma (|e|nGj = |f2B|/27rc is 
the Goldreich- Julian charge density). The magnetization pa- 
rameter characterizes the quotient of electro-magnetic flux 
to particle kinetic energy flux near the surface of an object. 
The so-called cr-problem, i.e., the problem of transformation 
of electro-magnetic energy into particle kinetic one, appears 
while one trying to explain an effective particle acceleration 
in the magnetic field. 

Indeed, the flow in the vicinity of different objects is 
assumed to be strongly magnetized at its origin. In spite 
of the lack of observational data near the surface of ra- 
dio pulsars, theoretical modeling predicts that the wind 
in this region has a composition cr ^ 1 (Michel 1 991; 
iBeskin. Gurevich fc Istominlll993) . The same can be said 
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about the AGNs iBeeelman. Blandford fc ReejlTosll . As 
for the Lorentz factor 7, at least for the blazars, the defi- 
ciency of a soft X-ray bump in their spectrum, which would 
be produced by the Comptonized direct radiation from the 
disc, is a reason to exclude the pre sence of particles with 
7 > 5 in the near zone of an object JSikoraet aJl200ED . On 
the other hand, at large distance from a pulsar, observations 
and modeling allow us to determine the magnetization pa- 
rameter (T ~ 10^'^ (Kennel fc Corornti 1984). Observations 
of quasars and act ive galactic nuclei give a ^ 1 and a <^ 1 
(|sikoraera]|200i). 

Up to now the axisymmetric stationary MHD approach 
gave the inefficient p article acceleration beyond the fast 
magnetosonic surface fBeskin, Kuznet sova fc Rafikovlll998l : 
pBoaovalov 1997; Bcaclman fc Li 1994). This seemed to be 
a general conclusion for any structure of a flow, however a 
lack of acceleration in the super sonic region was rather the 
consequence of monopole fleld l|Spitkovskv fc Aronj|2004l : 
iTh omson. Chang, fc Qua taer1^'2004^. The reconnection in 
a striped wind was discussed as one of possible accel- 
eration scenarios JCoronitillT990L iLvubarskv fc Kirkll200ll : 
iKirk fc Lvubarskvll200lF 7 however, this mechanism can ac- 
count only for the non-axisym metric part of a Poynting 
flux. Besides, for the Crab pulsar. ILvubarskv fc Kir'^ (1200 
had shown that effective acceleration may take place only 
beyond the standing shock. Another acceleration process 
can be connected with possible restriction of the longitu- 
dinal current, and thus the appearance of a light surface 
|E| = |B| at the flnite distance from a central object. In 
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this case the effective energy conversion and the current clo- 
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In this work we are trying to solve the problem of 
particle acceleration self-consistently within the approach 
of stationary axisymmetric magnetohydrodynamics. We re- 
gard particle inertia to be a small disturbance to the 
force-free flow. This allows us to linearize stream equa- 
tion and to find the disturbance to the magnetic flux, 
which corresponds to the finite mass of particles. Given 
that we can find the growth of a Lo r entz f actor. In the 
work of iBeskin. Kuznetsova fc Rafikovl l)l998h the Michel's 
monopole solution was taken as the zero approximation to 
the fiow. For this structure of the magnetic field, the Lorentz 
factor was found to be a^^^ on the fast magnetosonic surface, 
which was located at the finite distance unlike the force- 
free limit. Beyond that singular surface the acceleration 
turne d out to be ineffec tive. Treat ing the prob l em nu mer- 
ically, iBoeovalovl 1^^) and later iKomissarovl ||2004) also 
has got inefficient acceleration and coUimation for monopole 
outflow. Now we took the fiow near another force- free solu- 
tion, i.e., Blandford's paraboloidal magnetic field jBlanfordl 
Il97fill . The obvious difference of this zero-approximation in 
comparison with monopole solution is a well-coUimated fiow 
even in the force- free limit. 

In Section 2 we remind the trans-field and Bernoulli 
equations describing the stationary axisymmetric ideal out- 
fiow. After formulating the problem in Section 3, it is shown 
in Section 4 that the fast magnetosonic surface (FMS) is 
located at the distance rp — Ri^{g /OY^^ from the central 
object, the Lorentz factor 7 being (crS)^'^^ on it. This result 
is consistent with the standard value 7 a^^"'' . In Section 5 
it is shown that the disturbance of magnetic flux due to the 
finite mass of particles remains small in the subsonic region, 
the Lorentz factor growing linearly with a distance from the 
rotational axis within the FMS. In Section 6 we show that 
the supersonic fiow in the paraboloidal geometry is in fact 
one-dimensional, and it is possible to treat the problem nu- 
merically. In the supersonic region the growth of Lorentz 
factor remains linear with the distance vj from the rota- 
tional axis, reaching the value a at the distance z ~ a^Rh 
from the equatorial plane. This corresponds to almost the 
full conversion of the Poynting fiux into the particle kinetic 
energy flux. Finally in Section 7 we discuss some astrophys- 
ical applications. 



2 BASIC EQUATIONS 

Let us consider a stationary axisymmetric MHD fiow of cold 
plasma in a flat space. Within this approach magnetic field 
is expressed by 

B^^-^e,. (2) 

Here ^'(r, 0) is the stream function, I(r,9) is the total 
electric current inside magnetic tube '^(r,6) = const, and 
vj = r sin 6 is the distance from the rotational axis. In this 
paper we put c = 1. Owing to the condition of zero longi- 
tudinal electric field, one can write down the electric field 
as 



V*, 



(3) 



where the angular velocity f2p is constant on the magnetic 
surfaces: 57f = Q.y(^)- The frozen-in condition E-|-vxB = 
gives us 

u=-'B + 'yQ.Y-aeui, (4) 
n 

where u is four- velocity of a flow, and n is the concentration 
in the comoving reference frame. Function 77 is the ratio of 
particle fiux to the magnetic field flux. Using the continuity 
equation V(nu) = 0, one gets that 77 is constant on magnetic 
surfaces as well: r) — ri{'i/). 

Two extra integrals of motion are the energy flux, con- 
served due to stationarity. 



(5) 



and the z-component of the angular momentum, conserved 
due to axial symmetry, 

-^(*) = 77" + M'^^ip- (6) 
2n 

Here fi is the relativistic enthalpy, which is a constant for the 
cold flow. The flfth integral of motion is the entropy s(>I'), 
which is equal to zero for the cold flow under consideration. 

If the flux function 'i' and the integrals of motion are 
given, all other physical parameters of the flow can be de- 
termincd using the following algebraic relations dCamenzindl 
1^986; Beskin 1997): 

I _ L-Qftu^E 
2^^ ~ 



1 
1 



X2' 



1 {E-nFL)flFzu'^ -M'^L 



where the Alfvenic Mach number M is 



(7) 
(8) 
(9) 

(10) 



To determine , one should use the definition of Lorentz 
factor 7^ — = 1 which gives the Bernoulli equation in the 
form 

K _ 1 >I^(V*)2 



■co'^A^ 647r4 
Here 



2 2 
+ ^l r] . 



(11) 



K = uj'^{E- Q.FLf{l - ^i-co'^ - 2M'^) + M'^ivj^E^ - i^). 

The cold transonic flow is characterized by two singu- 
lar surfaces: the Alfvenic surface and the fast magnetosonic 
surface (FMS). The flrst is determined by the condition of 
nulling the denominator A in the relations (I7J-(|3J. FMS 
can be defined as the singularity of Mach number's gradi- 
ent. Writing equation IIH in the form 



647r'* K 



647r 222 



and taking the gradient of it, we can get 
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Here 
D = 



A 1 



= (vW — + 2^W^''^' 

and the operator V^. acts on all quantities except for M'^. 
The regularity conditions 

D ^0; Nr = 0; Ne ^0 

define the position of FMS and the relation between the 
integrals of motion on it. 

Finally, t he stream equ ation on the function ^'(r, 6) can 
be written as jBeskinlll997l) 



(1 



+ 



where 



(12) 



(13) 



Operator d/d'i/ acts only on the integrals of motion. The 
stream equation 1121 contains the magnetic flux function >!' 
and four integrals of motion; £■(*), L{'^), r?(*), and fiF(*), 
i.e., it has the Grad-Shafranov form. 




Figure 1. Blandford's force-free paraboloidal solution. Here the 
thick arrows represent the volume currents. The shadowed region 
is the working volume Hp const. The dashed line represents 
the light cylinder. 



3 THE PROBLEM 

Our goal is to determine the characteristics of a flow in the 
paraboloidal magnetic field. For this reason it is convenient 
to use the following orthogonal coordinates: 

X = r{l- cosO); Y = r{l + cos6); if. 

Here X stands for the certain magnetic surface in the force- 
free Blandford's paraboloidal solution, Y is the distance 
along it, and is the azimuthal angle. The latter does not 
appear in the equations because of the axial symmetry of 
the problem. The fiat metric in this coordinates is 



5xx = 



X + Y ^ 
iX ' 



gvY 



X + Y ^ 



g^pxp — XY. 



The n, Blandford's for ce-free solution can be written 
down as (iLee fc Park^l2004^ : 



" ^i + n|(x)f5 



/o(*) 



C»f(X)X 



(14) 



(15) 



where fip is the arbitrary function of and C is a constant. 
In particular, for fip = const we have 



S2p \ 



(16) 



We assume that in the vicinity of a central object the 
particle energy flux is much smaller than that of electro- 
magnetic fleld. In this case it is possible to consider the 
contribution of particle inertia as a small disturbance to the 



quantities of the force-free flow. Thus, in the first approx- 
imation we can get from (1121 the linear equation on the 
disturbance and solve the problem self-consistently. 

As was already stressed, for cold plasma the problem 
is characterized by two singular surfaces: Alfvenic and fast 
magnetosonic ones. Consequently, we need to specify four 
boundary conditions on the disc surface D llBeskia.l997i) . 
For simplicity, we consider the case 

nF(*)iD = const = fip, 

7jD = const = 7in < cr^''^, 

'?('I')|d = const = 77, 



*Id = 



^In f fipcc 4- J(yiYVoY + 1 

Sip v 



The condition f2p = const naturally restricts the region of a 
fiow under consideration (see Fig. 1). Since we assume that 
the magnetic field is frozen in the disc, we must consider 
the fiow only when f2pX|i5 = f2pcc7 < 1. In fact, we should 
use the inequality f2pX <^ 1. Then, Michel's magnetization 
parameter a for our problem can be defined as 

_ Sa _ Cflp 
[IT] Snurj 

Here Ea is a kind of energy amplitude: 
C^pX C^pX 



(17) 



E 



= Ea ■ 2npX. 



(18) 



Under our assumptions cr 2> 1, and we can introduce the 
small quantity e = . Besides, we will be mainly interested 
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in the flow far from the light cylinder, and so we have another 
limitation Q^XY > 1. 

We shall seek the stream function of the problem in the 
form 



<f{X,Y)=^o{X)+ef{X,Y), 



(19) 



where ef{X, Y) is the disturbance to find. The function of 
the angular momentum L may in general be different from 
the function Lq. For L the following expression may be writ- 
ten: 

,aL(*o) 



L(\[') =L(*o)+e/- 



9*0 



L(*o)+e^. 



(20) 



4 FAST MAGNETOSONIC SURFACE 

In order to find the position of FMS one can rewrite the 
Bernoulli equation IIH in the form 



4 , r, 3 

q +2q 



2 2 

urn 



1 



g +2g 



2 2 



(21) 



= 0. 



S2 ' f^pt^jS 

Here by definition q = , 



and e' = E{'^) — QfL(vI') = A"77(l ^ f^Fron^) ^ 7in/iJ? 
const. It is easy to check that for the force-free solution 

1 - cose 



(22) 



< 1 



(23) 



2 4 

for the small angle 9, i.e., in the whole region near the ro- 
tational axis. Remember that ^ = for Michel force-free 
monopole outfiow. 

To show that the quantity q is much smaller than unity, 
one can write © assuming a 3> 1 and ^^tiP ^ 1. Solving 
this equation for q, we get 

.= ^-^. (24) 

Thus, q is approximately equal to the ratio of particle kinetic 
energy to the full energy of the fiow, i.e., q <^ 1 for the 
magnetically dominated flow. 

As a result, one can rewrite 12111 in the form 

{e'f 



2g - U + 



1 



2 2 



£2 Q.Y-VJ'^ 



0, 



(25) 



where the terms and q were omitted due to their small- 
ness. 

Fast magnetosonic surface corresponds to the in- 
tersection of two roots of equation l)25|l . or to the 
condition of discriminant Q bein g equal to zero 
jBeskin. Kuzn etsova fc Rafikovl Il998h . The regularity 
conditions give dQ/dr = and dQ/d9 = 0, or dQ/dX = 
and dQ/dY = 0. For equation 1251 discriminant Q is 
expressed by 

i44 -ii22/ - \3 

1 Tj 1 1 n Tj 



16 £4 



16 27 E'^ 



The condition of root intersection Q 
as 



(26) 

can be rewritten 




Figure 2. Position of the Alfvenic surface (long dashed line) and 
FMS (short dashed lines) for paraboloidal magnetic field. The 
self-similar flow can be realized only outside the working volume. 



= 3 



and the first regularity condition dQ/dY as 
di 1 



Y 



dY 



Q.IXY' 



(27) 



(28) 



Taking approximately ^ « Y • d^/dY we get for the position 
of FMS 



'^F(e) ^ i?L ( J) 



1/2 



(29) 



where _Rl ~ f2p^ is the radius of a light cylinder. 

Values of q{rF,9) and 7f = ^{rp,9) due to condition 
Q = do not depend on the sum (^ + and on the 

FMS are equal to 



q{rF,9) = 
7F = {a9) 



erf 

1/2 



(30) 
(31) 



Again we confirm that q <^ 1 since a9 = Q,pXY\r-p ^ 1. 
These results are valid when 



. E 



EQ,fzu 



(32) 



i.e., in the region where electro-magnetic energy is greater 
than kinetic one. On the FMS this region is defined by the 
angle 9 changing from 7inCr~^ to a^^^. Here the greatest 
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value of 8 is given by the condition i}pX — 1, which cor- 
responds to the boundary of the working volume. On the 
other hand, for 9 < 7^1°""^ one can get 



rp « _Rl - 



7F = 7ii 



(33) 
(34) 



As we see, along the rotational axis particle energy remains 
the same as near the origin. 

The position of the FMS on the rotational axis can be 
evaluated independently. Indeed, for — the condition 
D = coincides with the condition A = 0, i.e., the position 
of the fast magnetosonic surface coincides with the Alfvenic 
surface on the axis. Assuming that Up ~ ■^in and using the 
definitions Q a^nd 11711 one can obtain 



It gives 



rp ~rA 



-Rl, 



(35) 



(36) 



coinciding with 13311 . This distance is much larger than the 
appropriate radius rp « {cr /yi^f^^ Rl for monopole mag- 
netic field. But the transverse dimension of the working vol- 
ume 117 ~ a^^^Ri^ remains the same. We see as well that on 
the FMS within the working volume the value of Lorentz 
factor changes from 7in to a^^^ . The greatest value of 7 is 
equal to the corresp onding value for the monopole structure 
iBeskin. Kuznetsov a & Rafikov 1998). 

Finally, as one can see in Fig. 2, the shape of the FMS 
extends along the rotational axis. Thus, the working volume 
cannot be described within the the self-similar approach as 
in this case the singular surfaces are to have conical shape 
9 = c onst iBtanford fc Pavn&.1982; .Li. Chiueh fc Beeelmanl 
Il992l). 



5 SUBSONIC FLOW 

For the inner region of the fiow we shall write down the 
stream equation with the small disturbances to the functions 
^ , L, and E. We shall treat the quantities q and ef as being 
of the same order of smallness. In the zero approximation 
one can get the equation 



d 



X + Y dX 



X 



dx 



+ 



+■ 



XY d^o/dX 



aL(*o) 
dX 



(37) 



= 0. 



Clearly, it has a solution II6II for I/(5'o) = -Lo(^'o)- On the 
other hand, in the first approximation we have 



X + Y 



X + Y 



1 647r* 



'dX 



X 



9*0 
dX 



X 



dq a*o 



dX dX 



X 



df d^f df 

dX dY^ dY 



+ 



2 XY^ d-Ho/dX dX 



+i 647r-* , df/dX dL^-^o) ^ 



(38) 



XY [d^o/dxy dx 



647r'' 



9(L(*o)0 



XY d^a/dX dX 



+ 



+- 



647r** 



= 0. 



XY ilp d^o/dX dX 

The integral L depends on the variable y as * = + £,/ 
depends on it. 

Before we proceed to solve equation 1381 we can eval- 
uate the ratio ef/'^o on the fast magnetosonic surface. In 
order to do this, we need to express ^ as a function of ef: 



^ ttC \X dXj 



(39) 



Given the order of ^ on the FMS by 12811 . one can get 

ef 1 



ttCX 
Thus, 



-5 « 1' 



(40) 



(41) 



where cr^ 3> 1 as we are interested in the flow structure 
outside the light cylinder. Hence, our disturbance proved to 
be small in comparison with the force-free solution up to the 
fast magnetosonic surface. Finally, from 1381 we can get an 
equation on the function ef: 

d^f 



+-kC 



(2qX 



X- 



dq_\ 

dXj 



+ 



(42) 



: 0. 



To obtain q, we shall make a natural assumption that 
q and ^ grow monotonically from correspondingly and 
near the origin of the fiow to {firj/Ef^^ and l/r2| zu'^ on 
the fast magnetosonic surface. Thus, in 12511 one can ne- 
glect the terms g'' and ^q^ in comparison with (firi/E)^ and 
q'^ /{flpi^)^ correspondingly. After that the solution of 1251 
is expressed by 

'^-2-Ax) • 

It is necessary to stress that in equation 1421 we can ne- 
glect the derivatives over Y in comparison with the deriva- 
tives over X. Here we take into account our assumptions 
that Y > X, so that df/dX ~ //X by the order of mag- 
nitude. Inside the working volume this means that we can 
neglect the curvature of field lines in our problem. Indeed, 
one can find that 



2e 



(xy)i/2 



a/ 9V 
dY 9y2 



(44) 
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formally writing the expression for curvature k for the im- 
plicit function '^o{X{x,y)) + £f{X{x,y),Y{x,y)) = const. 
Here the term {X^^^ /Y^^^) corresponds to the inverse cur- 
vature radius of the force-free magnetic surfaces. Thus, the 
curvature term does not play any role in the force balance 
on the magnetic surfaces in the paraboloidal magnetic field, 
which was quite different for the monopole magnetic field. 
In the latter case the curvatu re term played the leading role 
in the asymptotic region OBeskin fc Okamota , 2000 'l. 

After substitution of the found function q{X, Y) into 
14211 ■ one can find for the disturbance of the stream function 

ef ^—npiXY^^^ (45) 
(J 

with the Lorentz factor being equal to 
7 = QpiXYf^^ = nprsmd. 



(46) 



Thus, in the subsonic region the Lorentz factor grows lin- 
early with the distance from the axis and reaches on the fast 
magnetosonic surface the value (crS)'^/^, which corresponds 
to the result found in the previous Section. In this sense one 
can say that the solution in the inner region of the flow can 
be continued to the fast magnetosonic surface. 



6 SUPERSONIC FLOW 

In order to solve the problem in the supersonic region we 
need to emphasize some features of the paraboloidal config- 
uration of magnetic field. 

(i) The character of the flow may change in the vicinity 
of the singular surface. 

(ii) It was shown above that for paraboloidal magnetic 
field the curvature of field lines does not play any role in 
the force balance inside the working volume flpX < 1. This 
allows us to consider the fiow as one-dimensional. 

(iii) Positions of the fast magnetosonic surfaces in the 
paraboloidal field and in the cylindrical field coincide. 

Let us clarify the last point. The position of the fast magne- 
tosonic surface for the magnetized cylindrical jet immersed 
agnetic fi 

jBeskin fc MalvshkiDir2000l) 

Bext 



(47) 



where to ~ is the distance from the rotational axis, and 
'i'jct is the total fiux inside the jet. On the other hand, ac- 
cording to our definition of a I17II for paraboloidal flow 



E = 2firjanFX 



47r2 



(48) 



As a result, we shall deflne the magnetic flux inside the 
region fipX < 1 as ^Pjet, and thus C is expressed by 



C = 



(49) 



In this case the expressions for a for two flows coincide. 
Taking poloidal paraboloidal fleld as the external field for 
the one-dimensional flow 
C 



Br, = 



2z' 

one can get the position of the FMS 



(50) 



(5) 



1/2 



(51) 



which coincides with (1291 . 

Hence, the flow becomes actually ID in the vicin- 
ity of the FMS, not to say about the supersonic region. 
For this reason we can c onsider th e supersonic flow as 
one-dimensional. But unlike ISeskin fc M alvslikin (2000) we 
would use the paraboloidal magnetic field -Bp(z) II50I I out- 
side the working volume as an external one. It gives us the 
slow z-dependence of all the values. 

For the cylindrical flow the integrals of motion near 
the axis are the same as the integrals of motion in our 
paraboloidal problem: 



r2F(\l/) = const, 



(52) 
(53) 



77(5') = const, (54) 
£(*) = 7inW + ^fL + flpL, (55) 
where e' = const. Introducing non-dimensional variables 



(56) 
(57) 

one can rewrite equations 11H - 112II as a set of ordinary 
differ ential equations for y and M l|Besk in fc MalvshkirJ 
l2000h: 



(1 



M 



^'11) 



+x^{r,,, + 2yf 



(7in + X -1) 



Ay'^M^ 



dx 



2 2 

^ Tina: 
x^ 

— (1 - 



2xM 



{1-x^ 



2M^)+ 



(58) 



-flxM^ 



(1- 



-+ 



(59) 



2;3(1 



We want to emphasiz e that in the work of 
iBeskin fc MalvshkirJ i2000l) these differential equations 
were applicable only for the inner part of the jet, where the 
integrals of motion could be expressed as the linear functions 
of 5' and Qf = const. In fact, in the case <^ x^ (and 
that is in the range of parameters we are interested in) the 
solution for Bp is the same for an arbitrary function Sip (5' 
the e quation I58I I can be rewritten as (.Beskin fc Malvshki 
12001 



B^zu) = 



(60) 



As E{'^) is proportional to the Q,p{'^), the particular expres- 
sion for the angular velocity is not contained in the equation 
for Bp. This ensures the continuity of the solution even in 
the region where the condition f2p — const does not hold. 

Analytically, from the set of equations (I581 - (I591 one 
can get the following results: 



(i) X < 7in 

M'^ = MI^ const, y='^' 



2Mr 



(61) 
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i.e., the poloidal magnetic field is approximately constant, 
(ii) X > 7in 



(a) > jI 

..2 Mo 2 ,2; 
J\4 = — —X , y (X m — . 

7in 



7ii 



(62) 



>r2 



i.e., the poloidal magnetic field decreases as B „ oc ^ 
Chiueh, Li & Bcgolman .1991 : lEichler. ilQQSi : iBogovalovl 



(b) Ml « 7f„ 
M = X, y = 



7in 2 

a; , 



(63) 



7in 27^2 

i.e., again Bp ~ const. 

Using the connection 
Lorentz factor 

Iirj X 

Then, for x 3> 7i n and <C the foll owing linear de- 



lO' 10' 10" lo' 10° 

Distance from the disk, z 



10 



10 



'jjj.ri/E, one gets for the Figure 3. The border of the working volume xjct as a function of 

z for 7in = 8 and a = 10^. Its form remains almost paraboloidal. 
All the distances are given in the units of the light cylinder. 



(64) 



pendence is valid jBeskin fc Malvshkinlliobal : 



7 ^ 



(65) 



Let us find the distance along the axis until which the linear 
growth of the Lorentz factor continues. In order to do this 
one should write the constant poloidal magnetic field which 
defines M^: 



*on| dy *of^l7in 



27ro-x dx 2TvaJ\4o 



(66) 



As this magnetic field should be equal to the outer one, we 
get 



Cr7in-RL, 



(67) 



and the greatest Lorentz factor near the boundary of the 
working volume is 

7 = ('^7in)''". (68) 

For z > afinRh we shall perform numerical calculation 
in the intermediate region between two limits A^o ^ 7iii 
and A4o ^ 7in. The boundary condition is the equality of 
the magnetic fiux and of magnetic field. Thus, for every z 
we need to find the proper value of Mq that would allow 
us to make internal poloidal magnetic field on the border of 
the working volume be equal to the external paraboloidal 
magnetic filed, as the internal fiux becomes equal to the 
total flux of a jet ^'jot- In this case the border of the working 
volume Xjct is defined during the numerical integration. We 
expect it to be 



Xjct{z) 



1/2 



(69) 



The integration shows that it remains almost paraboloidal 
(see Fig. 3). Such an integration gives the linear growth of 
Lorentz factor until about a (see Fig. 4 and Appendix). 

Thus, the Lorentz factor 7 grows linearly with the dis- 
tance from the axis reaching the value a near the border 
of the working volume for z = a^Ri^. This corresponds to 
the transformation of about a half of electro-magnetic en- 
ergy into the kinetic one. The maximal Lorentz factor in this 
problem is 2a: 



10 



10= 



O 4 

5 10 

N 

c 

S 10 



10 
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10^ 










X 10^ 









































10 



10 



10 



10^ 



10= 



10" 



Distance from the axis, Xj^j 



Figure 4. The growth of the Lorentz factor 7 near the border 
of the working volume as a function of the distance from the 
rotational axis (which is given in the units of the light cylinder) 
for different values of cr (10^, 10*, 10^) in a logarithmic scale. The 
linear growth of Lorentz factor continues until ■y ^ cr. 



E_ 



= Jin + 2(7, 



(70) 



which follows from the definition ||HJ. Thus, in the 
paraboloidal magnetic field the effective particle accelera- 
tion can be realized. 

One can see that the condition for the particle acceler- 
ation from the work l)Vlahakisll200l 

holds for our geometry. Indeed, in our case this condition 
written for the border of the working volume 'I' = ^/jot trans- 
forms into the 



_9_ 



< 0. 



(72) 
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As the function x'j^^{z)/z decreases along the field line (see 
Fig. 3), the condition ^ is fulfilled. 

Here we should point out the main difference of our 
problem from the work in which the Michel's monopole so- 
lution was taken as the first approximation. 

(i) Even the force-free fiow is well coUimated already. 

(ii) The curvature term of the paraboloidal problem does 
not play any role in the force balance, which allowed us to 
treat the flow as one-dimensional in the supersonic domain. 

(iii) The considered fiow is supported from outside by the 
much greater fiux than the one which is contained in the in- 
ner region. This allows the flow to widen and, consequently, 
to accelerate outflowing particles inside the working volume. 



lO' 10' 10" lo' 10° 

Distance from the disk, z 



10 



10 



7 MORE REALISTIC MODELS 

In the previous section we were discussing a model in which 
the rotational velocity fip of the field lines, which were as- 
sumed to penetrate a horizon of a black hole, in the working 
volume was constant, and it was equal to zero outside the 
working volume. This allowed us to regard the vacuum field 
outside the working volume as confining the fiow. In this 
case the particle are accelerated effectively up to a. However, 
such model does not seem to be realistic since we assume the 
non-rotating disk. 

In this section we shall discuss two models. They would 
have three distinct regions. The first one is the working vol- 
ume with the given law for rotational velocity, where we 
assume the MHD fiow of the electron-positron plasma pro- 
duced by the Blandford-Znajek process. The second one is 
the region of the, presumably, slow ion wind originating from 
the disk rotating with Keplerian velocity 



np(y) = QF(cr) 



3/2 



(73) 



where we choose the parameters so that the function SlF(y) 
would be continuous on the border of the working volume. 
The third region is the vacuum field where the disk angular 
velocity drops almost to zero. For our convenience we put 
Qf being equal to zero at the distance where the Keplerian 
velocity drops ten times less than the rotational velocity 
near the axis. So we still have some portion of the vacuum 
magnetic fiux to support our fiow configuration. 

The first model would be still characterized by fip = 
const inside the working volume. For the second one we em- 
ploy the following rotational velocity Qf, operating in the 
working volume: 



(l-a;2)[l-hln(l-ha;) 



{41n2 -I- 1 - 3;2 + [1 - x2 - 2(1 + a;)]ln(l + x)} ' 



(74) 



Here x = cos&h, and On is the spherical angle at the hori- 
zo n labeling the field line. It was first found in the work 
of lBlanford fc Znaiekl I^T^) for the paraboloidal magnetic 
field threading the slowly rotating black hole. So the second 
model must be the most realistic. 

For these two models we would perform the numerical 
calculation for the following set of the ordinal differential 
equations: 



Figure 5. The border of the working volume Xj^t as a function 
of z for = 8 and a = 10^ . This curve is the same for both the 
first and the second models. Its form is only slightly compressed 
in comparison with parabola. All the distances are given in the 
units of the light cylinder. 



\dxj 



2 2 



{A-M^)+ 



J7|(0)' 



^2(0) 



(75) 



7iU 



"l(y) 
ni.(o) 



x' -\ 



dx ^ Ax^ ^f{0) 



xM^ Qliy) 



A 



i-yl - 2A) + 



M dy ( X 



dnliy) 



(76) 



n|(o) 



dy 



The boundary condition again is the equality of the mag- 
netic flux and of magnetic field at the border between the 
vacuum fleld and the internal flow with the variable rota- 
tional velocity fip. We shall be interested in the character- 
istics of the fiow in the working volume, assuming that the 
wind from the disk is governed by the MHD equations. 

For the first model we see the same law for the particle 
Lorentz factor as a function of the x on the border of the 
working volume as the one, we have gotten in the previous 
section. But now the border is compressed greater in com- 
parison with its shape in the previous section (see Fig. 5), 
so the growth of the Lorentz factor along the axis is slightly 
slower. 

The first question with the second model is whether 
the fiow in subsonic region is greatly affected by the non- 
constant angular velocity of the field lines. It is easy to show, 
that the linearized equation 1421 for the given function fip 
can be written as 



^2 d'f df 



£L ( 

dx \ 



d'^f df 



+-kC 



2qX (l + 



-(l-iM)=0, 



(77) 



where AI 



dnl 



is much smaller than the unity for 



QfX < 1. 

The numerical integration shows again that the border 
of the working volume is compressed greater than it was for 
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Figure 6. The growth of the Lorentz factor 7 near the border 
of the working volume as a function of the distance from the 
rotational axis (which is given in the units of the light cylinder) 
for different values of cr (10^, 10"', 10^) in a logarithmic scale 
is presented for the second model. Again the linear growth of 
Lorentz factor takes place, but it grows until approximately 0.5(t. 



the vacuum field just outside the working volume. Besides, 
the growth of the Lorentz factor, although sustaining the 
shape of the previous curve ln7(lna;), stops at smaller values 
than in the case of a constant rotational velocity (see Fig. 6). 
this results from the diminishing of the Qf near the border 
of the working volume: 



7n 



= 2(7 



n|(o)- 



(78) 



As the n|(cr)/n|(0) 0.28 for the Blandford-Znajek solu- 
tion, so does the Lorentz factor at the border of the working 
volume. We may conclude that the decrease of the function 
Qf generally suppresses the effectiveness of particle acceler- 
ation. 



8 DISCUSSION AND ASTROPHYSICAL 
APPLICATIONS 

In reality, there are at any way two reasons why the ideal pic- 
ture under consideration may be destroyed. First of all, for 
large enough 7 the dra g force Jdrag 7'^ can be i mport ant. 
As was demonstrated in lBeskin. Zakamska fc Soil i2003l) . for 
monopole magnetic field it takes place when the compact- 
ness parameter Ia becomes larger than a^^^. Second, the 
particle acceleration continues only until the paraboloidal 
poloidal magnetic field can be realized. This can be stopped 
by the external magnetic field, say, by the galactic chaotic 
magnetic field -Bgai ~ 10~^ G. The fiow inside the border 
QpX = const widens, while the external paraboloidal mag- 
netic field -Bext is greater than the Bgai . When they equate, 
the expansion stops and the particle acceleration ceases. The 
distance Xst fi"om the axis at which _Bgai = Bext is given by 
— Bo/Bc^f Here Bo = ^'jetf2|. is the magnetic field 
strength in the vicinity of the central object. Thus, 



7max = Min{xst, cr}. 
The relation 
Bo N 



7 ; 



/_Bo_ 

V -Boxt 



(79) 



(80) 



was already obtained in lBeskin fc MalvshkinI ()2000l) . As can 
be easily seen from 1791 . whether the Poynting flux will be 
transformed into the particle kinetic energy fiux depends on 
the value of a and Bcxt- 

We now consider several astrophysical applications. 



8.1 Active Galactic Nuclei 

For AGNs the central engine is assumed to be a rotating 
black hole with mass M ~ 10^ Mq, R ~ 10" cm, the total 
luminosity L ~ 10''^ erg s"\ Bo ~ lO" G. In the Michel 
magnetization parameter a Q 



io'^x-^i.bJ—) 



V c / 



(81) 



the main uncertainty comes from the multiplication param- 
eter A, i.e., in the particle number density n. Indeed, for 
an electron-positron outflow this value depends on the ef- 
ficiency of pair creation in the magnetosphere of a black 
hole, which is still undetermined. In particular, this process 
depends on the density and energies of the photons in the 
immediate vicinity of the black hole. As a result, if the hard- 
photon density is not hig h, then the multiplication parame - 
ter is smafl (A ~ 10- 100:lBeskin. Istomin fc Parie\^ Jl992l) : 
iHirotani fc Okamotd Jl99^ ). In this case for (QR/c) ~ 0.1- 
0.01 we have a ~ 10^ - 10^^. Knowing a we can estimate 
the maximal Lorentz factor 7max 17911 . In the presence of 
the external magnetic field of typical value Bext = 10~^G, 
7max = 10"" <C CT, so Only the small part of the Poynting 
fiux can be transformed into the particle fiux. On the other 
hand, if the density of photons with energies > IMeV is 
high enough, direct particle creation 7 -I- 7 — > -I- e~ re - 
sults in an increase of the particle density JSvenssoJliii) . 
This gives a ~ 10 - 10^ In this case the Lorentz factor is 
7max = o", and the energy tra nsformation can be efficient. 

R ecently the work of iGracia. Tsinganos. Bogovalovl 
12005*) has been published in the electronic arXiv with an 
idea close to the one presented here (the central relativis- 
tic flow with the disk wind surrounding it), although with 
different formalization. In that work the comparison of the 
observational data for the jet opening angle with the model 
opening angle is presented. Here we can repeat the same 
for our paraboloidal inner fl ow. The distance fo r the an- 
gle range from the work of iBiretta et al.l ll2002l) lies well 
inside a sphere with the radius z = 7?L7inO", so we shall use 
the strictly paraboloidal form of the border of the work- 
ing volume. We have one free parameter in this case: the 
value of f2po which we chose to be equal to 8.2 • 10~*. The 
result is presented in the Fig. 7. Although the analytical 
curve fails to explain the points at an angle « 33 degrees, 
it fits well the outer region and (surprisingly) the innermost 
point. The Lorentz factor for such distances is « 7in, so the 
jet possesses the same moderate Lorentz factor f avour ed in 
the works oflBiretta et all Jl995l) . lBiretta et al.l Jl99gl) and 
ICramphorn et al J J2004I) . 
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Figure 7. The observational data of the opening angle for M87 
and the border of the worlting volume for Op(0) = 8.2 ■ 10^** s^^ . 



8.2 Radio Pulsars 

For radio pulsars the central engine is a rotating neutron 
star with M Mq, R ~ l(f cm, and Bo ~ 10^^ G. In 
this case the magnetization parameter a ~ 10^-10®, corre- 
sponding to relativistic electron-positron plasma, is known 
with rather high accuracy (see, e.g., Bogovalov ( 199^). The 
Lorentz factor in the presence of an external magnetic field 
is 7max = 10* ~ cr. 

This estimate is valid, of course, only for the aligned ro- 
tator. In a case of an orthogonal rotator the situation with a- 
problem is more explicit, as in this case the Goldriech- Julian 
charge density must be the (QR/c)^^'^ smaller near the polar 
caps than that for the axisymmetric magnetosphere. Thus, 
it is natural to expect that the longitudinal current flowing 
alon g the open field lines would be p roportionally smaller 
too iBeskin. Gurevich fc Istominlll993t) . But, consequently, 
toroidal magnetic field will be smaller than the poloidal elec- 
tric field in the vicinity of the light cylinder. On the other 
hand, it is known for the Michel's monopole solution that 
in order to remove the light surface to infinity, the toroidal 
magnetic field must be of the same order as the poloidal 
electric field on the light cylinder. If the longitudinal current 
j does not exceed by {D,R/c)~^^'^ times the quantity PqIjC, 
where Pq'j is an average charge density on the polar cap 
with X ~ 90° (for the typical pulsars this factor approach 
the value of 10^), the light surface for the orthogonal rota- 
tor must be located in the vicinity of the light cylinder. In 
this case the effective energy conversion and the current clo- 



iBeskin. Gurevich & Istomin 


I1993I: IChiueh. Li & BeeelmanI 


ll998l:lBeskin & Rafikovll200C 


)■ 



field of So ~ 10^^ G (which is necessary to explain the to- 
tal energy release) the magnetization parameter a is small 
(cr < 1 — 100), because within this model the magnetic field 
itself is secondary and its energy density cannot exceed the 
plasma energy density. Would it be not so, i.e., the magnetic 
field would be prior to the particle energy fiow, we could for- 
mally apply our estimate. Having the external magnetic field 
to be the order of 10^^ G, we can get the standard value for 
the Lorentz factor 7max ~ 10^. However, in this case it is 
hard to explain of the value of a being the order of 10^. 



9 CONCLUSION 

We have gotten the characteristics of the flow in the 
paraboloidal magnetic field within the approach of station- 
ary axisymmetric magnetohydrodynamics. To simplify the 
problem, we assumed that in the strongly magnetized flow 
with cr ^ 1 a particle inertia could be described as a small 
disturbance to the force-free flow. As the zero a pproxima- 
tion t he solution with paraboloidal magnetic field iBlanfordI 
Jl976h was taken. 

The position of the fast magnetosonic surface is found 
to be rp ~ Ri^, with the Lorentz factor changing 

from 7in to a^^^ on it. The disturbance e/ to the stream 
function * is e/ = nC{nlXYy^^/a < *o inside the FMS. 
As to the Lorentz factor 7, it grows linearly with the distance 
from the axis. 

On the fast magnetosonic surface the structure of the 
flow may change significantly. It is implicitly confirmed by 
the fact that the characteristics of the flow, which we got 
under the assumption of the small disturbance in the super- 
sonic region, is not in the agreement with the results on the 
FMS. However, in our problem the curvature term does not 
play role in the force balance on the magnetic surface, and 
the positions of the FMS in the cylindrical and paraboloidal 
flows coincide. These facts allowed us to regard the problem 
as one-dimensional and to perform numerical calculations. 

As a result, we got the further growth 7 = {z/R^)^^^ 
of the Lorentz factor until it reaches the value of a near the 
border of the working volume for z ~ a^Rh from the equa- 
torial plane. This corresponds to almost the full conversion 
of the Foynting energy flux into the particle kinetic one. 

For the more realistic model with three regions — the 
jet with the Blandford-Znajek rotational velocity, the disk 
wind, and the vacuum field — the particle acceleration is 
suppressed, but still the maximal Lorentz factor is of the 
order of a. 

We want to emphasize th at our solution cannot be de- 
scrib ed in a self-similar way JContopoulojll995l : IVlahakisI 
I2OO4I) as we assumed the constant angular velocity of the 
magnetic surfaces near the rotational axis. This structure 
cannot be considered within the self-similar approach. 



8.3 Cosmological Gamma-Ray Bursts 

For cosmological gamma-ray bursts the central engine is rep- 
resented by the merger of very rapidly orbiting neutron stars 
or black holes with M ~ Mq, R ~ 1 0^ cm, and tota l lu- 
minosity L ~ 10^^ erg s~^ (see, e.g., iLee at all ||200(]|) for 
detail). On the other hand, even for a superstrong magnetic 
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APPENDIX A: THE ANALYTICAL SOLUTION 
FOR THE ONE-DIMENSIONAL FLOW 

Let us consider the set of equations 1581 - 15911 . For x 2> 7in 
and y 2> 7in the system can be rewritten as 

2 



dx - ■■■ • x^{l — x^-M^) 
Making the following substitutions 



^ _ 1 



y 



(Al) 
(A2) 

(A3) 



and treating z as a new variable, one can get the following 
set: 



dq 
dz 
dt 
dz 



t{l 



2z 



(A4) 



(A5) 



It is convenient to transform this set into the one second- 
order differential equation 

d^q 1 + q dq 3 f dq^ 
dz^ 



= 0. 



2z dz 

Knowing its solution, one can get the function t{z) as 
tiz) = -^. 



(A6) 



The solution of IIA6II gives 
2a^/b{q-qi){q-q2) 



l/v/T^ 



y = 
7 = 

Here 



ab{q - qi){q - q2) 



12 



2 



qi){q - q2) 



9 + 1 



1 - y]~~b 1 + \/r~~6 

gi = r , 92 = - 



(A7) 



(A8) 



(A9) 



(AlO) 



(All) 



6 ' 6 ' 

a and b are constant, < a, —00 < 6 < 1. For b < the 
maximal value of q is equal to unity. On the other hand, 
for 6 > the value of q is not limited. Thus, we would 
consider only the case 6 > 0. For the Mq < 7hi th^ solution 
is represented by the lower branch, and the case M/q > 7j^ 
by the upper branch of the graph (see Fig. Al). 

Let us find the values of a and b. In order to do this 
we shall regard two limits of this solution: for TVfo «C 7in 



and lor Ml > 7in. 



The condition Mq ^ 7;^ holds when 



q <Si qi < g2 , so we get the known solution 

- 2 1 /g2^ ^''^ 
y = Ax = -p — 
4a \ qi 



(A12) 
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2a(b)/ 



ab 



y(q),x(q) 



depending on the quotient Mq /'yf^, which changes from 1 to 
(compare this result to the numerical calculation in Fig. 6). 
Lorentz factor is expressed in this case by (see Fig. 7) 



7 oc a; 



(A20) 



So, if we assume that the effective grow of the Lorentz-factor 
continues until the exponent in llA20t is equal to 1/2, we 

can estimate the maximal Lorentz-factor by the value — 

^ 2 

(compare it with the Fig. 4). 



This paper has been typeset from a Tp^X/ ffl^jX file prepared 
by the author. 



Figure Al. The functions x(q) and y(q). 



— — — 'la\/hq 



(A13) 



Equations A = 7in/(2Mo) and B — M^jy^^ are not inde- 
pendent and give the following connection between a and 
6: 



Ml / 1 + vr~6 



27ii, 



(A14) 



When Mq 2> 7?n <i ~ const = Mq /7?n, which corresponds to 
the asymptotical approach of g to g2. Thus, for g2 = Mq /7in, 
one can get 



6 = 



{Mlhlf ■ 



(A15) 



Thus. 



Ml ( Ml 



(A16) 



Finally, for Mq 7?„ (see the lower branch in the 
Fig. Al) the expressions 1631 hold. In the intermediate re- 
gion Ml « 7?„ (the upper branch in the Fig. Al) noth- 
ing can be said except the general formulas <A8IA9ll . For 
Ml ^ 7hi (the upper branch in the Fig. Al, q 52) we 
can get the refined expressions for y{x) and M'^{x). Decom- 
posing the functions lASII - IIAQII near 52 as g — 52 = S, we 
get 



y 



Thus, 

^2.2(VT=6-l)/(yT=6-2) 



y 



Ax' 



(A17) 
(A18) 



(A19) 



where A is constant. Thus, for q = const the expression for y 
is not Ina;, but the exponential function of x with the index 



